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ABSTRACT 

We test whether the model proposed by Esin et al. to explain X-ray observations of the black 
hole soft X-ray transient GRS 1124—68 in outburst can also explain the optical lightcurves of a 
similar object A0620— 00. We show that to reproduce the observed X-ray to optical flux ratio in 
A0620— 00, we need to assume X-ray irradiation of the outer disk that is significantly in excess 
of what is expected for a standard planar disk. With enhanced irradiation, the Esin et al. 
model can reproduce the optical evolution of A0620— 00 in outburst very well. Though we find 
that optical observations of GRS 1124—68 also imply enhanced X-ray heating of the outer disk, 
the irradiation appears to be a factor of ~ 3 weaker in this system than in A0620— 00. This is 
surprising, since GRS 1124—68 has a larger disk. We speculate that enhanced irradiation may 
be due to disk warping, and that the degree of warping differs between the two binaries. 

Subject headings: accretion, accretion disks - black hole physics - plasmas - radiation 
mechanisms: thermal 



1. INTRODUCTION 

A number of known X-ray sources exhibit dramatic 
outbursts in which their bolometric luminosity in- 
creases by many orders of magnitude. A fraction of 
these systems, termed soft X-ray transients (SXTs), 
is believed to contain binaries in which a more mas- 
sive compact star is accreting mass from a less mas- 
sive late-type companion (see Tanaka &: Shibazaki] 
1996 and van Paradijs & McClintock 1995 for the 
most recent comprehensive reviews, though several 
new SXTs were discovered in the last few years). This 
subset of X-ray transients is especially interesting be- 
cause the low mass of their secondary companions has 
in several cases allowed us to measure the mass func- 
tions of the compact primary to be above 3Mq, thus 
providing the first (and so far only) unambiguous ev- 
idence for the existence of stellar mass black holes 
(see iMcCUntock 1998| and [Charles 1999| for the most 



The behavior of black hole SXTs (BHSXT) during 
their outbursts shows a number of common features. 
A "typical" lightcurve consists of a fast rise (r ~ 
few days) followed by a longer (r ~ month) exponen- 
tial decline back to the system's quiescent luminos- 
ity, sometimes with one or more secondary maxima, 
though not all BHSXTs follow this standard prescrip- 



tion (see Chen, Shrader, & Livio 1997 for an excellent 
inventory of observed lightcurves). The X-ray spec- 
tral evolution of these systems during outbursts is 
less well studied. There are many indications that 
the majority of sources start out fairly soft near the 
peak (hence their name), i.e. with the total emis- 
sion dominated by a component with T ~ few keV, 
and soften further during the decline, before becom- 
ing suddenly very hard (the total energy output is 
dominated by ~ 100 keV photons) several months af- 
ter the peak (e.g. p?anaka &: Lewin 1995|; [Tanaka 



up - to - date list of black hole SXTs). 



fc Shibazaki 19^ ). However, there are several ex- 



ceptions; e.g. GS 2023-^338 and GRO J0422-F32 re- 



^Chandra Fellow 
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mained hard throughout their outburst, though the 
latter displayed a very typical exponential lightcurve. 
Very little is known about the X-ray spectra of BH- 
SXTs in quiescence, due mainly to their faintness; 
only V404 Cyg has a well determined spectrum (see 
Narayan, Barret fc McClintock 19971) . 

It is now generally accepted that the outbursts in 
BHSXTs are caused by a thermal instability in the 
accretion disk, the same mechanism that is believed 
to cause outbursts in dwarf novae (see e.g. Cannizzc| 
1993|-] |Lasota 1996| ; |Can"nizzo 1998| ). This model relies 
on the fact that in an optically thick accretion disk 
with a central temperature Tc ~ 10^ K (the ioniza- 
tion temperature of hydrogen) the opacity is a very 
strong function of Tc. Both cold neutral disks and 
hot fully-ionized disks are stable, but when the cen- 
tral disk temperature is close to the critical value, 
thermal instability forces the disk to undergo a limit 
cycle. In quiescence the disk is cold and its mass ac- 
cretion rate is considerably below the mass transfer 
rate from the companion. In this way the disk accu- 
mulates mass until increasing density and tempera- 
ture force it into a hot ionized state, characterized by 
a very high (often on the order of MEdd) rnass accre- 
tion rate, thereby starting an outburst. Of course this 
high mass accretion rate can only be sustained until 
the surface density in the disk is depleted sufficiently 
to force a transition back into the cold neutral state, 
after which the system settles back into quiescence. 

A detailed numerical model of such an outburst in 
a BHSXT is currently under const ruction (see e.g. 



ing the decline from the outburst peak to quiescence, 
Esin, McClintock, fc Narayan] (1997, hereafter EMN) 



Dubus, Lasota &: Hameury 1999a| ; [Hameury et al 
1998|]|. It is already clear however, that illumination 
of the outer disk by X-rays produced near the black 
hole plays a very i mportant role (e.g. van Paradijs &j 
McClntock 1994 Ivan Paradi.is 199^; [Lasota 1996|). 



King Sz Ritter (1998) in a recent paper plausibly 
argued that such irradiation would naturally pro- 
duce a prototypical exponential lightcurve. These 
authors also attribute the departures from the stan- 
dard lightcurve shape to the size of the disk, claiming 
that in larger systems irradiation is unable to heat up 
the outer disk to temperatures above 10^ K. In such 
a case a cooling wave setting off from the outer disk 
edge would shut off the mass flow early on, producing 
a short outburst with a linear decline. This sugges- 
tion seems to be supported by observations, as shown 



by Shahbaz, Charles fc King (1998), as well as b y de- 
tailed theoretical modeling ( Dubus et al. 1999a|) . 

Though the disk instability model outlined above 
offers a plausible explanation of the general lightcurve 
shapes seen in outbursts of BHSXTs, the origins 
of the spectral evolution of these systems are still 
very much debated in the literature. In an at- 
tempt to unify different "spectral states" seen dur- 



have proposed a model for the prototypical BHSXT 
GRS 1124-68 (also Nova Muscae 1991). In their sce- 
nario, luminous soft spectra (with characteristic pho- 
ton index T ^ 2.5) seen near the peak of the out- 
burst, when the system is in the so called Very High 
and High states (hereafter VHS and HS respectively), 
are modeled by a cool optically thick accretion disk 
( ^hakura &: Sunyaev 1973D topped with a hot o pti- 



cally thin corona (e.g. Haardt & Maraschi 1991). In 
the less luminous Low state (hereafter LS), when the 
observed spectra are very hard (with T ~ 2.0 — 1.5), 
as well as in the Quiescent state, the accretion flow 
assumes a different configuration: a hot and opti- 
cally thin inner flow (a sort of diskless corona) is 
surrounded by an optically thick outer disk with an 
inner edge anywhere between several tens and tens 
of thousands Schwarzschild radii from the accreting 
black hole. 

EMN model the hot inner part as an optically thin 
advection-dominated accretion flow (ADAF) discov- 
ered by Ichimaru (1977) , which has been extensively 
discussed in the literature in the last six year s (see 



Narayan, Mahadevan, fc Quataert 1998| and [Kato 



Fukue & Mineshige 1995 for latest reviews). An 



ADAF is a hot and stable accretion flow solution al- 
ternative to the standard Shakura & Sunyaev (1973) 
thin disk. It is characterized by electron tempera- 
tures above 10^ K and optical depth ^ 1, which are 
ideal conditions for producing the hard spectra seen 
in the low state through thermal Comptonization of 
the seed synchrotron and thin disk photons. Since an 
ADAF solution only exists for mass accretion rates 
below some critical value. Merit; it suggests a natu- 
ral explanation for the transition from the hot flow -|- 
outer thin disk configuration at low M to thin disk 
-|- corona at M > Merit- 

A model involving an advection-dominated flow at 
accretion rates well below Merit surrounded by an 
outer thin disk had already been shown to explain 
well the puzzling properties of quiescent BHSXTs 



(e.g. 


Narayan, McChntock, & Yi 1996 


|al. 1997 
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Hameury et al. 1997). In their 



demonstrated that the scenario outlined above works 
for systems in outburst; it reproduces very well both 
soft and hard X-ray lightcurves of GRS 1124—68 as 
well as the details of its spectral evolution. However, 
though GRS 1124—68 had excellent X-ray coverage. 



the optical (e.g. 


Bailyn 1992 




King, Harrison, McNa- 


mara 1996|; pclla Valle, Masetti & Bianchini 1998) 



vations, which are the only diagnostic of the outer 
disk properties, were limited to the first 150 days of 
the outburst. At later times EMN could obtain use- 
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ful ccnstraints only on the innermost regions of the Kolb &: Burderi 1996). Finally, the size of the accre- 



accretion flow, the properties of which are well con- 
strained by the observed X-ray emission. Particularly 
unfortunate is the fact that the important High-Low 
state transition occurred during a period with almost 
no optical coverage. According to the EMN model, 
during this transition, the source of the X-ray emis- 
sion shifted from the flat disk to a quasi-spherical 
corona, causing a considerable change in the angular 
distribution of radiation incident onto the outer disk. 
This effect could potentially be seen in the optical 
light curve. 

Though we cannot remedy the lack of optical cov- 
erage for GRS 1124—68, in this paper we attempt to 
do the next best thing, i.e. we attempt to comple- 



ment the work of EMN by using the extensive multi- 
wavelength photometry available for the outburst of 
another BHSXT, A0620— 00, to infer the properties 
of the outer disk in these two systems. The qual- 
ity of the X-ray data available for A0620— 00 (though 
superior for its time; see Kuulkers 199j ) renders it 
unusable for detailed modeling similar to that done 
for GRS 1124—68. However, in sections and §|3|we 
will demonstrate that the two systems are sufficiently 
alike in their binary parameters and outburst proper- 
ties to justify a direct use of the model proposed for 
GRS 1124—68, with only minor parameter changes. 

In §^ and §|5| we describe the modeling of the optical 
data and present our results, which are then further 
discussed in We conclude with a summary in §^. 

2. A0620-00 AND GRS 1124-68 AS BINARY SYSTEMS 

The remarkable similarity between GRS 1124—68 
and A 0620— 00 in their orbital parameters (Bailyn 



19921 ; [Remihard, McClintock, k Bailyn 19921), ap 



pearance of the secondary (Remillard et al. 1992) 



and of the accretion disk in quiescence ( Orosz et al 
1994), has been well documented in the literature. 



To emphasize how closely these two systems resem- 
ble each other, we show in Table |^ a summary of 
their latest binary parameters, gathered from the lit- 
erature. 

The properties of BHSXT outburst lightcurves de- 
pend largely on the mass of the primary, which deter- 
mines the overall energy scale for the system, mass 
transfer rate from the companion, and the size of the 
accretion disk. The latter is very important since 
it probably dictates the value of the mass accretion 
rate at the peak of the outburst as well as the over- 
all hghtcurve shape (e.g. [King fc Ritter 19"98t) . It 
is clear from Table ||, that the data are consistent 
with the black holes in GRS 1124-68 and A0620-00 
having similar masses. Since their orbital periods are 
short, this argues that the mass transfer rate can- 
not be very different in the two systems (e.g. Kingj 



tion disk, iiout, can be estimated from other orbital 
parameters, assuming that it fills a fixed fraction of 
the primary Roche Lobe (we take it to be 80%). We 
compute i?out using the paczyiiski (1971)] formula: 



Rout = 2.82 X 10i°P'/if 

X (a ?!8 — n 91r.o-n~ 



+ (2-1) 

(0.38 -0.2 log g) cm. 



where Porb is the binary orbital period measured 
in hours, q = M^/M is the mass ratio, and m = 
M/Mq is the mass of the black hole in solar units. 
In calculating accretion disk size for A0620— 00 we 
adopted m = 4.5 (see ^ and q = 0.067, while 
for GRS 1124-68 we took m = 6 and g = 0.133. 
The resulting values of i?out for both objects with 
the corresponding uncertainties are shown in Table 
|l|. Since the orbital period of A0620— 00 is shorter 
than that of GRS 1124—68, the accretion disk of the 
former is somewhat smaller in physical units, though 
in Schwarzschild units the two disks have roughly the 
same sizes. 

3. X-RAY OUTBURSTS OF GRS 1124-68 & A0620-00 

Given that these two binaries resemble each other 
so closely, it is not surprising that their properties 
in outburst are similar as well, as we demonstrate 
below. In Figure |^ (upper panels) we have plotted 
the observed soft X-ray lightcurves of A0620— 00 and 
GRS 1124—68. Both systems show a fast rise and 
a characteristic exponential decline after the peak, 
interrupted by a small secondary and a somewhat 
larger tertiary maximum. (Though it is not clear 
from this data that A0620— 00 reaches a local maxi- 
mum near day 200, other instruments showed that the 
X-ray flux indeed declined shortly after that point, 
see Kuulkers 1998|.) The rates of the exponential de- 



cline after the primary maximum in the two objects 
are very close; for A0620— 00 we have r ~ 32 days, 
and for GRS 1124-68, r ~ 35 days. 

It is considerably more difficult to compare spectral 
evolutions of the two BHSXTs, since all published 
spectra of A0620— 00 were taken during the first three 
months after the onset of the outburst. These obser- 
vations all show fairly soft spectra with characteristic 



peak temperatures around a few keV (e.g. Carpen- 



ter et al. 1976; Citterio et al. 1976; Poxsey et al 
1976 ; Long &: Kestenbaum 1978| ) confirming that this 



system indeed looked like GRS 1124-68 in the VHS 
and HS. After about ~ 100 days, we have to rely 
on the broadband data, so in Figure || (lower panels) 
we have plotted ratios of count rates in the hard and 
soft X-ray bands for both objects. The figure shows a 
nearly identical slow decrease of the hardness ratios 
for both systems, followed by an abrupt hardening 
just prior to the rise to the tertiary maximum. The 
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latter feature suggests that a High-Low state transi- 
tion is taking place. 
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Fig. 1. — (a) Soft X-ray lightcurve and hardness ratio 
of A0620-00 obtained with SAS-3 CSLA (1.5 - 6keV) and 
CSLB (6 - 15keV) detectors (courtesy K. Plaks, J. Woo 
and G. Clark; adopted from Kuulkers 1998). (b) Soft X-ray 
Hghtcurve and hardness ratio of GRS 1124—68 obtained with 
Ginga LAC detectors (courtesy of K. Ebisawa). Note that 
error bars are not plotted in either panel. 



So far in our comparison of the X-ray lightcurves 
and hardness ratios we used raw count rates, since 
this is the information most readily available for 
A0620— 00. Though in general a questionable pro- 
cedure, we believe it is justified here because Ginga 
LA C and SAS-3 contained detectors with very similar 



responses in the energy bands of interest ( Buff et al 
19771; [Turner et al. 1989| ). We conclude that the ex- 



isting data strongly supports our hypothesis that the 
X-ray outbursts of these two BHSXTs showed almost 
identical evolution. However, our goal is to adapt 
the EMN model to A0620— 00, and since we cannot 
convert the existing A0620— 00 data to physical flux 
units (or at least not enough of it to be of use) for 
comparison with the model, we have to obtain some 
sort of quantitative mapping between GRS 1124—68 
and A0620— 00 lightcurves. If this can be done, then 
we can assume that the model for GRS 1124—68 
(with properly adjusted parameters, as deduced from 
the lightcurve mapping) reproduces the behavior of 
A0620— 00 in the X-ray band and proceed to test it 
against the optical data. 

HjIVHN" adopted the following mapping relation be- 
tween the time after the outburst maximum (T, mea- 
sured in days), mass accretion rate (M) and the tran- 
sition radius between the inner hot flow (ADAF) and 
the outer thin disk (Rtr)'- 



M, 



peak 



M 



+ TRtr In 



Rtr 

3i?fr 



(3-1) 



of Mpeak and the transition radius is at its minimum 
value of 3Rg (the marginally stable orbit for a non- 
rotating black hole), i.e. most of the mass is accreting 
via a thin disk. As T increases, the mass accretion 
rate decreases exponentially (with a characteristic e- 
folding ti me scale Tj^-j) as the di sk empties out, as 



argued by [King &: Ritter (1998) . In this regime, a 



dominant soft component originates in the thin disk 
and hard X-ray emission comes from a hot corona 
above the disk. According to the relative importance 
of the hard X-ray component, |EMN identified the re- 
sulting spectra with the VHS (stronger high energy 
emission) or HS (very weak high energy emission). 

Once M drops below the critical value for the for- 
mation of an ADAF, Merit, the disk begins to evap- 
orate starting from its inner edge, forming the hot 
accretion flow and thereby increasing the transition 
radius (on an e- folding time scale tr^J. This is ac- 
companied by a dramatic change in the X-ray spec- 
trum. Since now most of the emission originates in 
the very hot, rarefied gas filling the region inside Rtr, 
the spectrum quickly becomes very hard, with most of 
the emission coming out in 100 keV photons, which 
are produced via inverse Compton scattering of seed 
synchrotron and disk photons. This regime, when M 
through the inner region stays near its critical value 
(fueled by the evaporation of the disk, s ee §|6[) and 
Rtr increases with time, was identified by EMN with 



the Intermediate state (hereafter IS). In this state, 
spectral and temporal properties of BHSXTs shift 
from those characteristic of the HS (soft spectrum) 
to those of the LS (hard spectrum). Finally, at some 
point (perhaps when Rtr reaches its value in quies- 
cence) mass accretion rate through the inner region 
begins to drop again and the system moves from the 
LS into quiescence. 
In general, there is no reason to expect the param- 



eters in equation (3-1) to be exactly identical for two 
systems, even as otherwise similar as those consid- 
ered in this paper. However, in the scenario pro- 



At T = 0, mass accretion rate is at its maximum value 



posed by EMN, there is one quantity which must 
remain relatively constant from system to system, 
and that is rhcrit = -Merit /MEddi where MEdd = 
1.39 X lO^^m gs~^ is the Eddington mass accretion 
rate for a black hole of mass m = M/Mq. Therefore, 
we expect that in Eddington units the X-ray lumi- 
nosity of BHSXTs in the IS is more or less constant. 

Based on this fundamental assumption, in Figure 
1^ we have plotted soft X-ray lightcurves for our two 
BHSXTs, where the A0620— 00 flux was scaled to ob- 
tain the best match in the section of the lightcurves 
corresponding to the intermediate state, as marked on 
the figure. The fluxes are normalized to Crab units 
to remove the effect of different instrumental effec- 
tive areas. Then the scaling factor r] = Fa/Fgrs = 8 
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must be due only to the difference in the black hole 
masses, binary inclinations and distances to the two 
systems. Roughly, we have 

_ Fa _ / ruA \ f cosiA \ f dcRS \ ^ 
^ Fgrs \mGRs) \cos icRs) Wa / 

(3-2) 

Note that though the dependence of the observed flux 
on cos i used here is strictly valid only for an optically 
thick disk, the disk emission dominates soft X-ray 
flux during early part of the IS. Moreover, near rhcrit 
the emission from ADAF also escapes preferentially 
along the polar axes, where the gas density is lower 
(see Narayan fc Yi 1995 and Narayan et al. 19971) . 
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Fig. 2.— Soft X-ray lightcurves of A0620-00 {SAS-3 
CSLA; 1.5 - 6keV) and GRS 1124-68 {Ginga LAC 1.2 - 
9.3 keV) plotted in Crab units. Different spectral states seen 
during the outburst of GRS 1124—68 are indicated in the fig- 
ure. The count rate for A0620— 00 has been multiplied by a 
factor of 1/8, to match the curves in the IS, when, according to 
the EMN model, mass accretion rate is on the order of rhcrit. 
Note the close correspondence between the lightcurves. They 
differ only in the time of the secondary maximum and in the 
width of the tertiary peak. 



' EMN* adopted mcRS = 6, igrs = 60° and dcRs = 
5kpc and showed that their model reproduces the 
X-ray data of GRS 1124—68. Given these values 
for GRS 1124—68 and the constraints on the dis- 
tance to A0620-00 (see Table |), equation (|^ fa- 
vors a somewhat low black hole mass for A0620— 00. 
For niA = 4.5, we obtain dA = 1.4 kpc or 1.5 kpc, 
depending on the value of the f{M). Larger val- 
ues of niA require placing A0620— 00 at an even 
greater distance, apparently inconsistent with exist- 
ing observations. An even smaller black hole mass 



in A0620— 00, e.g. niA = 3.5, is formally consistent 
with all the data; however, there is a preference to- 
wards a somewhat higher value. Therefore, in our 
model for A0620— 00 we adopt the largest black hole 
mass that is still consistent with the distance to the 
system: uia = 4.5; ia = 65°; dA = 1.4 kpc. 

The values of the three parameters for A0620— 00 
introduced in equation (S^l) can be read off Fig- 
ure |2[ By a direct comparison of the two outburst 
time scales, we obtain Tj^,j = 32 days, quite close to 
Tj^.j = 35 days for GRS 1124—68. Assuming that at 
day 200 Rtr = 10^'^i?g (determined from observations 



of line in quiescence, see [Narayan et al. 1996 and 



references therein), where Rg = 2.95 x 10 m cm is 
the Schwarzschild radius, the length of the interme- 
diate state implies tr^^ = 10 days (as compared with 
TR^^ = 8 days for GRS 1124-68). Finally, the ra- 
tio of the peak luminosity to that in the IS is the 
same in both systems, which implies that the peak 
mass accretion rate in A0620— 00 was the same as in 
GRS 1124-68, i.e. rhpeak = Mp^^^/M^dd = 3. Inci- 
dentally, similarity of rhpeak in these two systems sug- 
gests that like GRS 1124-68, A0620-00 was in the 
VHS near the peak of the outburst. The available 
spectral information, though limited, also supports 
this conclusion. 

One must keep in mind however, that our assump- 
tion of constant rhcrit has been demonstrated to be 
valid (in theory) only for accreting black holes with 
zero spin; and it is possible that the value of rhcrit de- 
pends on a/M. If the black holes in A0620— 00 and 
GRS 1124—68 have very different spins, the proce- 
dure of matching the two lightcurves described above 
may not be valid. On the other hand, Zhang, Cu: 



Chen (1997)1 and |Cui, Zhang fc Chen (1998)| argue 



that GRS 1124—68 does not contain a rapidly spin- 
ning black hole, and since A0620— 00 shows similar 
X-ray spectra, it is probably not likely to have much 
spin either. 

4. MODELING UBV LIGHTCURVES OF A0620-00 



The model of EMN] underpredicted the optical flux 
of GRS 1124-68 in the VHS and HS by a factor ~ 3, 
though outer disk irradiation was self-consistently 
included in their calculations. This is not surpris- 
ing, since it has recently been demonstrated in sev- 



eral papers (e.g. see Dubus et al. 1999b and refer- 
ences therein) that the standard Shakura fc Sunyaev 
(1973)1 thin disks do not provide enough flaring at 



the outer edge to reproduce the observed levels of ir- 
radiation. In fact, self-consistent disk models show 
that the outer parts of the disk (even when irradia- 
tion is taken into account in the disk height calcu- 



lations) are convex, rather than concave ( Tuchman 



[Mineshige, fc Wheeler 1990| ; Pubus et al. 1999b|) , i.e 
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they are shielded from the X-rays emitted by the in- 
ner disk. Since in the luminous spectral states (VHS 
and HS) most of the X-ray emission originates in the 
thin disk, rather than the corona, this places very se- 



vere constraints on how important irradiation can be 
in planar disks. 

As we will demonstrate below, the ratio of the opti- 
cal to X-ray fluxes in A0620— 00 is even greater than 
that in GRS 1124—68. Clearly then, in order to re- 
prodi jico tho oboorvod optical omiDDion from tho for 



mer s yatcm, wc need to aooumc that the ourfacc of the 
outer disk is considerably more curved than predicted 
by the standard model. An alternative explanation, 
involving an extended hot corona which can illumi- 
nate even a concave outer disk, cannot apply when 
the X-ray spectrum is dominated by the disk emis- 
sion and fails to explain the statistics of disk eclipses 
in low-mass X-ray binaries ( |de Jong, van Paradijs <^ 
Augu lsteijn 1996| ) . 

The well-known equation describing irradiating 
flux at radius R due to a compact central source of 
X-rays is given by (e.g. phakura fc Sunyaev 1973 ; 
van Paradi.is 1996| ; [King fc Ritter 19981) : 



nLx(l 



47ri?2 



R 



d\nH 
d\nR 



1 



(4-1) 



where Lx is the total central X-ray luminosity, a is 
the disk albedo, and H is the height of the disk pho- 
tosphere above the orbital plane. The power index 



in |roduced by King fc Ritter (1998)|, is equal to 



n, 

1 when the irradiating source radiates isotropically 
(e.g. as in the case of a corona); and n is set to 2 
when the outer disk is irradiated by the inner disk, 
which emits primarily in the direction perpendicular 
to the disk surface. The factor of n in the numerator 
(significant only in the case when n = 2) is neces- 
sary to ensure the correct normalization of the total 
emission from the thin disk. 
We know Lx from analyzing the X-ray data (see 
, which also allows us to distinguish whether n = 1 
or n = 2 (or a combination of both) is applicable. 
Unfortunately, the other two parameters, H{R) and 
a(R), are not well known. The former because, as 
we argued above, the standard answer cannot re- 
produce the observations. The latter because pub- 
lished theoretical calculations of X-ray irradiation, 
from which a value of a can be derived for a given 
incident spectrum, generally average over all incident 
angle s (e.g. Magdziarz fc Zdziarski 1995|; Zycki et al 



1994|-] [White, Lightman fc Zdziarski 1988[) , while here 



only grazing incidence X-rays are relevant. In fact, 
the values of a inferred from observations of persis- 
tent low-mass X-ray binaries with neutron star pri- 
maries ( de Jong et al. 199(: ; Kallman, Raymond &| 
Vrtikk 1991) are very close to unity, quite different 



from the generally accepted range, a = 0.1 — 0.2 

We simplify the problem by lumping the two un- 
certainties together as a single irradiation parameter 
C{R) ( Pubus et al. 1999b| ; [Esin, Lasota fc Hynei ' 
19991 ), so that 



C{R) 



AttR^' 



(4-2) 



(Note that our C is defined through Lx, while Dubus 
et al. (1999b) define their parameter through Mc^ , so 
that the two differ by the factor r/ = 0.1, which rep- 
resents the efficiency of converting mass into X-rays 
for the thin disk.) Moreover, we assume that the only 
radial dependence of C comes from H(R), for which 
we adopt a rather ad hoc prescription: H oc i?^/^. 
In systems in outburst, the optical band generally 
falls in the Rayleigh-Jeans regime of the outermost 
disk radius. Therefore, the total optical emission is 
well described simply by the value of H at the outer 
disk edge, H{Ro^t), and the exact form of H{R) 
affects only the details of the optical spectra. We 
picked the prescription above derived for isothermal 
disks ( Vrtilek et al. 1990| ), solely to make the con- 
nection with previous literature on disk irradiation, 
though we know that the outer disks in A0620— 00 
and GRS 1124—68 are not isothermal. For the X-ray 
albedo we adopt (simply for lack of anything more 
appropriate) the angle averaged prescription for the 
X-ray irradiation of a cold unionized medium ( [White 
et al. 1988| ) which gives a 0.1 — 0.2, depending 



on the spectral distribution of incident photons. We 
believe that this is not a bad approximation, since 
grazing incidence irradiation is not likely to increase 
the value of the albedo above ~ 0.4, unless the outer 
disk has an ionized atmosphere with optical depth of 
order unity (D. Psaltis, private communication; see 
also Basko, Sunyaev fc Titarchuk 19"?^) . Thus the 
error in the quantity (1 — a) in Eq. ( [4-1] ) is at most 
50%, and not an order of magnitude, as implied by 
the albedo values derived by [de Jong et al. (1996) 



Nevertheless, because of the uncertainty in this quan- 
tity, the values for //(i?out) quoted below cannot be 
taken at face value; they must simply be interpreted 
as a measure of C. 

With the prescription for the treatment of outer 
disk irradiation outlined above, we use the lEMN 



model for X-ray emission to compute a series of spec- 
tra following the evolution of A0620— 00 during the 
decline phase of its outburst. Throughout the cal- 
culation the shape of the outer thin disk was held 
fixed, and the height of the disk at the outer edge, 
H{Rout)/RovLt, was adjusted to fit the overall nor- 
malization of the observed optical emission. The ir- 
radiating flux originating in the ADAF was computed 
using the ray tracing procedure described in detail in 
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EMN , which is roughly equivalent to using Eg. ( 4- 1 
with n = 1. We calculated the irradiating flux emit- 
ted in the inner parts of the thin disk using Eq. ( |4-1| ) 
with n = 2. To describe the emission from the thin 
disk we used the standard modified black-body spec- 
trum. The effective temperature at each radius was 
computed by equating the black-body emission to the 
sum of the internal energy dissipation rate in the disk 
and the heating rate due to irradiation (as described 
in EMN). Note that in the outer disk, which is pri- 
marily responsible for the observed optical emission, 
X-ray irradiation of the disk surface completely dom- 
inates over the viscous heating. 

Optical Lightcurves of A0620-00 
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Model Parameters: 

BH = 4.5 
i = 65° 
D = 1.4 kpc 




V band (X = 5500 A) 
B band (A = 4400A) 
U band (A = 3650A) 




SO TOO rao 2tro 25t 
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Fig. 3. — (a) The observed specific luminosity of A0620-00 
in the three optical bands (Webbink 1978; Kuulkers 1998) is 
plotted as solid symbols together with the best fitting model 
lightcurves, obtained with H{Rout) ~ 0.12i?out. Note that 
most of the optical flux was generated by the reprocessing 
in the outer disk of the X-rays emitted by the inner disk 
and the corona (in the VHS and HS) or by the ADAF (in 
the IS and LS). (b) The variation of the mass accretion rate 
rh = M / M-Edd and the transition radius rtr = Rti/Rg with 
time in the model. 

The resulting optical lightcurves in the standard 
UBV optical bands, corresponding to the best fit- 
ting value of H{Ro^t) = 0.12i?out are shown in Figure 
|3|(a). In the same figure we plotted the observed op- 
tical fluxes for A0620-00 (from [Webbink 1978| and 
Kuulkers 19981) , dereddened for Eb-v = 0.35 ( |Wu et 



al. 198^ ) and converted to luminosity values assum- 
ing the distance to A0620-00 of 1.4 kpc. 

As in modeling the X-ray data, we used the 
mapping given by equation (p]) to construct our 



lightcurves. Figure |^(b) shows how two main model 
parameters, m = M/MEdd and rtr = Rtr/Rg, change 
with time during the decline phase of the outburst. 
The transition radius changes only during the IS, 
when mass accretion rate stays constant. Before and 
after the IS, rh decreases exponentially with time. 
Here we take the value of rtr at the end of the IS 
to be the same as that suggested by observations in 
quiescence ( Narayan et al. 1996| ). There is no di- 
rect observational support for this assumption from 
the studies of X-ray spectra, since the shape of X-ray 
emission stays essentially fixed for any value of rtr 
above ~ 30 (see e.g. Esin et al. 1998). However, a 



comparison of mass accreted during the IS with the 
total mass stored in the thin disk (see §^) favors a 
large change in rtr. 

The general correspondence between the model 
lightcurves and the data in Figure ^(a) is excellent. 
Taking into account the quality of the [EMN fit to 
GRS 1124—68 X-ray data, one cannot expect to pro- 
duce optical lightcurves based on this model that 
would be in better agreement with observations. 

Two points are important to emphasize. First, 
we reproduce the slower decay of the optical inten- 
sity (ropt ~ 75 days) relative to the X-ray intensity 
(tx ~ 25 days), which is a typical feature of opti- 
cal lightcurves in BHSXTs. Though optical emission 
comes purely from X-ray irradiation, cooling of the 
outer disk, which accompanies a decline in the irra- 
diating flux, causes the peak of the emission to shift 
towards the optical band. This effect compensates 
for the total decline in the reprocessed emission and 
considerably increases the e-folding time of the ob- 



served optical flux (see also the discussion in King & 
Ritter 199^ ). 

The second feature of the observed lightcurves 
which is reproduced well in the model, is the upturn 
during the IS (between days 130 and 200). In our 
model this is due to the change in the character of 
the irradiating flux. During the early parts of the out- 
burst (VHS and HS) X-ray emission comes primarily 
from the accretion disk, so that the irradiating flux is 
proportional to the square of the disk height (see Eq. 

1). During this period the irradiation parameter 
in our calculations is C(i?out) — 0.004. On the other 
hand, in the LS, when all the X-rays come from the 
isotropically emitting ADAF, the irradiating flux is 
linear in H/R. At this point, we have C(i?out) — 0.03. 
During the IS, the value of C changes gradually be- 
tween these two limiting values, causing an upturn 
in the lightcurve. This feature may be an additional 
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confirmation of our assumption that during tliis part 
of the outburst the transition radius is increasing. 

Though the overall normalization of the model 
lightcurves is determined by our choice of i?(i?out), 
it is interesting that the model produces the correct 
spacing between the V, B, and U lightcurves, indi- 
cating that our optical band spectra must agree with 
the data fairly well. However, because the functional 
form of H{R) was chosen arbitrarily, this may not be 
a very significant result. 

5. REVISITING GRS 1124-68 



As mentioned above, the original EMN model fell 
short (by a factor of 3) of reproducing the optical 
emission from GRS 1124—68. However, with the im- 
proved treatment of irradiation we now have no trou- 
ble in fitting those data, as demonstrated in Figure 
The UV data points were taken from Shrader & 
Gonzalez- Riestra (1993). The V band data are col- 



lected from Bailyn (1992)|; [King et al. (1996)| ; pella 



Valle et al. (1998); and other observations tabulated 
by Cheng et al. (1992, see references therein). [Chenj 
also obtained an HST spectrum of the sys- 
tem which we used to estimate the observed flux at 
3100A for day 127. All the data points were dered- 
dened using the color index Eb-v = 0.30 quoted 



by ^hrader fc Gonzalez- Riestra (1993) . The model 



lightcurves were calculated using the parameters for 
GRS 1124-68 given in |EMN| . 



Optical Lightcurves of GRS 1124-68 
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Fig. a.— Optical and UV photometry for GRS 1124-68 
are plotted as solid points. Errorbars on the UV points are 
quoted by Shrader & Gonzalez-Riestra (1993), and those on 
the V band points are all set to ±0.1 mag. The curves show the 
model lightcurves of GRS 1124—68 calculated for Tj^j — 35 days 
and Tijtj — 8 days (see EMN). The best fit is obtained with 
H{R 

out ) = 0.07R out • 

The most interesting point about the fit in Figure 
is the comparatively small best fit value of the disk 



height at the outer edge, if(i?out) = 0.07i?out- Using 
the fact that most of the data were obtained during 
the period when X-ray emission is dominated by the 
disk, we estimate C(i?out) — 0.0014. This value is 
considerably smaller than that derived for A0620— 00 
in §^ It is easy to see the origin of this discrepancy. 
A comparison of Figures ^(a) and Q shows that the 
optical luminosities of the two BHSXTs were nearly 
identical. However A0620— 00 was a factor of 1.6 dim- 
mer in X-rays (see §^) , implying that irradiation must 
have been considerably more important in this sys- 
tem than in GRS 1124-68. 

6. DISCUSSION 

Our results show that in order to reproduce the 
optical hghtcurves of A0620-00 and GRS 1124-68 
efficient irradiation of the outer disk is required. The 
disk height at the outer edge required to fit the ob- 
servations is a factor of ~ 2 smaller than H/R ^0.2 
found by de Jong et al. (1996) for some persistent 



LMXBs. However, this comparison of H/R values is 
complicated by the fact that the albedo values used 
in our calculations were almost an order of magni- 



tude smaller than those derived by de Jong et al 



(1996)| . In addition, we assumed a differe nt angU' 



lar distribution for the irradiating flux then de Jong 
et al., A better way to check our results would be 



through comparing the effective irradiation parame- 
ter at the outer disk edge instead of disk height, since 
it is really the former that determines the importance 



of irradiation. The results of de Jong et al. (1996) 
imply C(i?out) ~ 0.002 — 0.004, consistent with what 
we find in the VHS and HS. On the other hand, our 
calculations imply that in the IS the heating of the 
outer disk is almost an order of magnitude stronger 
because of the isotropy of the irradiating X-rays. 

The most puzzling result of our investigation is 
the fact that a smaller disk in A0620— 00 is heated 
considerably more than its larger counterpart in 
GRS 1124—68. This conclusion is nearly model in- 
dependent, since it follows simply from comparing 
the X-ray to optical flux ratios in the two systems. 
Of course our main assumption that the mass accre- 
tion rate during the IS is the same for A0620— 00 
and GRS 1124—68 may be incorrect, since we do 
not know precisely what mechanism is driving the 
evaporation of the thin disk at this stage of the out- 
burst, or whether the two black holes have the same 
spin. However, even abandoning this fundamental 
constraint would not reconcile the results for the two 
systems. If we allow m in the IS to be smaller than 
merit in A0620— 00, X-ray data allows higher m and 
smaller i for this system (which incidentally are also 
preferred by observations in quiescence, see Table 
|ID. The size of the disk in A0620-00 would not be 
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strongly affected since it depends only weakly on m 
(see Eq. p-l| ), so most of the effect would come from 
the change in z. In a lower inclination system the 
ratio of optical emission to X-rays is greater, as is 
observed in A0620— 00, but only if X-ray emission is 
close to isotropic. Thus, this would have some effect 
in the late IS and LS, where X-ray emission originates 
in an ADAF, but will not help in the VHS and HS, 
were X-rays are emitted predominantly in the disk. 

Thus, our comparison of A0620— 00 and 
GRS 1124—68 clearly implies that the disk in the 



The EMN] model does not explain the origin of 
the secondary maximum seen near day 55 in the 
lightcurves of A0620— 00 and around day 70 for 
GRS 1124—68. Among several possible explanations 
for this feature that were proposed in the literature 
(see Cannizzo 1998 and references therein), two sce- 
narios invoke an enhancement in the central mass ac- 
cretion rate triggered by irradiation of a new mass 
sour ce. Augusteijn, Kuulkers &: Shaham] (1993, see 
also [Chen, Livio &: Gehrels 1993|) propose the sec- 
ondary star as the source of the extra mass; [King &| 



form y is for some reason more 'curved', though in Ritter (1998)| invoke the outer disk, initially left cold 



all other respects the two systems are nearly identical. 
This suggests that the shape of the disk in a BHSXT 
may not be entirely predictable for a given set of 
parameters, or in other words, that it might depend 
on the past history of the system. The most obvious 



example of such an effect is disk warping (e.g. MaL_ 
oney, [Begelman fc Pringle 1996| ). If the outer disks in 
BHSXTs (and LMXBs in general) are warped, this 
can explain both high levels of irradiation and vari- 
ations between otherwise similar sources, like those 
we have examined here. 

One constraint we can place on disk warping comes 
from expected luminosity fluctuations due to the 
long-term precession of the warp. For the simplest 
warp geometry (m = 1 mode) we can estimate the 
ratio of fluxes observed when the warp is facing the 
observer and when it is facing away from the ob- 
server as cos {i — j3)/ cos {i + (3), where (3 is the an- 
gle between the normal to the disk surface at the 
outer edge and the binary axis, and i is the binary 
inclination. For our prescription H oc i?^/^, we have 
(3 = tan~^ (dH/dR) ~ 9°; in general for small warp 
amplitude /? is of the same order of magnitude as 
H/R. Thus, for our choice of parameters, the optical 
emission observed from a warped disk in A0620— 00 
must vary by a factor of ~ 2 on the time scale of 
several months (see e.g. Maloney et al. 1996 for 



radiation-driven warping), while the emission from 
GRS 1124-68 should vary by about 50%. (Note 
that no eclipses of either the secondary or the cen- 
tral X-ray source are expected in either A0620— 00 of 
GRS 1124-68, since in both systems P + i <. tt/2.) 
This number of course must be treated as only a 
rough estimate, as it is very sensitive to the exact 
shape of the disk surface. However, it is interest- 
ing to note that both optical and X-ray lightcurves 
of A0620— 00 show ~7-8 day oscillations of roughly 
the right amplitude (though the period seems to be 



too short for precession) during the IS (e.g. Tsunemi 



Mats loka & Takagishi 1977; Matilsky et al. 1976; see 
also Figures |2| and |^a) , and a hint of the same effect 
can be seen in the optical lightcurve of GRS 1124—68 
(Figure ^ see also the discussion in Kuulkers 1998| ) . 



by the heating wave. Both mechanisms would cause 
an increase in the lightcurve decay time after the 
peak of the outburst, where is essentially of or- 
der of the viscous time scale in the outer disk. In 
A0620— 00, a system with a much smaller disk, this 
time scale should be smaller as well, consistent with 
what is observed. 

The plateau before the tertiary maximum in the 
X-ray lightcurves, which we identify with the end 
of the IS, was much more prominent in A0620— 00 
than in GRS 1124—68. In our model this means that 
the phase of constant M and increasing Rtr lasted 
roughly ^ 20 days longer in A0620— 00. An inter- 
esting question to consider is why M in the X-ray 
producing region stays constant during this period. 
One obvious possibility is that this results directly 
from disk evaporation. We can write down the equa- 
tion for the mass accretion rate through the central 
ADAF during the IS: 



Madaf = Mdisk + (Rtr) 2TTRtrR 



tr- 



(6-1) 



This simply states that the mass flow through the 
ADAF is equal to the mass flow through the disk 
plus the rate of mass evaporation from the inner part 
of the disk. In outburst, we expect Mf^sk to de- 
crease exponentially with time. Then one can show 
that if Madaf is kept constant and the surface den- 
sity profile in the outer disk is taken to be that of a 
steady-state disk, the transition radius has to increase 
roughly exponentially with time. Note that we can 
assume that Madaf is constant with R and is equal 
to the central M, since the accretion timescale in the 
ADAF region is many orders of magnitude shorter 
that the timescale on which the quantities on RHS 
of Eq. ( |6-1|) vary, i.e. the ADAF adjusts practically 
instantaneously to the changes in the mass inflow at 
the transition radius. 

Though a realistic functional form of Rtrit) is prob- 
ably more complicated than a simple exponential 
adopted in Eq. ( ^-l|) (since it depends on the de- 
tailed evolution of the surface density in the outer 
disk, which can be obtained only from numerical disk- 
instability models), the main results of our paper 
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are not very sensitive to this detail. For our model 
lightcurves to match optical data for A0620— 00, it is 
sufficient for i?tr to increase in some fashion from 3i?g 
to ~ 102 - lO'^i^g during the IS. 

If most of the mass accreted onto the black hole 
during the IS comes from disk evaporation, it is im- 
portant to compare the mass accreted during this pe- 
riod to the mass stored in the thin disk at the onset 
of the IS. The former is simply 



7. CONCLUSIONS 
We have tested whether the model proposed by 



Merit Atis, 



(6-2) 



where Atis is the duration of the IS. To estimate the 
mass stored in the disk we make use of the well known 
fact that during the outburst the disk adopts a quasi- 
steady density profile, so that 



Rq 



Mdisk = 27r/ J:iM,,it,R) RdR, 



(6-3) 



where T,( Merit , R) is given by the standard outer disk 
formula from Shakura fc Sunyaev (1973) , and Rq is 
the value of the transition radius at the end of the 
IS. For A0620-00 we get 

W^-^R, 



- 3 X lO^'^g and 
6 X 1023(i2q/103-9i?g)5/4(0.l/a)4/5g, where 
is the value of the transition radius in qui- 
escence (see Narayan et al. 1996| ). Because these 
numbers are very approximate, the discrepancy be- 
tween them is not serious. However, we can safely 
draw the conclusion that without any other source of 
mass, during the IS the transition radius must have 
increased close to its quiescent value. 

According to this interpretation, the outer disk con- 
tains very little mass at the end of the IS plateau. 
This is consistent with the sharp drop in the X-ray 
(see [Kuulkers 199^ and references therein) and op- 
tical (see Figure |3^) fluxes observed at that point. 
Moreover, the last set of optical data points in Fig- 
ure §(a), observed on day 232, shows that the B band 
flux was considerably lower than both V and U band 
fluxes (this was originally pointed out by Lyutyi 1976| 
and [Lyutyi Sz Shugarov 1979 ), implying that perhaps 
the disk was becoming optically thin at this point. 

Though the IS plateau was somewhat shorter in 
GRS 1124—68, the same calculation for this sys- 
tem gives us very similar values: Mace ~ 4 x 10^^ g 
and Mdisk ~ 1024(i?q/lo3-9iig)5/4(o.i/a)4/5g. These 

numbers are clearly at odds with the results of ZyckiJ 
Done ! fc Smith (1998)| which seem to indicate that 
in GRS 1124—68 the transition radius was still fairly 
small at the end of the IS. 



Esin, McClintock, &: Narayan (1997) to explain the 
X-ray behavior of GRS 1124-68 during its 1991 
outburst can also reproduce the observed optical 
lightcurves of another BHSXT, A0620— 00. Our com- 
parison of the X-ray data for GRS 1124—68 and 
A0620— 00 strongly suggested that these X-ray novae 
must have followed the same spectral evolution dur- 
ing their outbursts and therefore the scenario devel- 
oped for GRS 1124—68 could be directly used in mod- 
ehng A0620-00. We showed that the [EMN| model 
modified to include a flared outer disk, can reproduce 
the shape of the optical lightcurves of A0620— 00. As- 
suming that the shape of the irradiated outer disk re- 
mains fixed throughout the outburst, this result pro- 
vides important evidence in support of our basic sce- 
nario of X-ray production, since it relies on the fact 
that the primary source of X-rays changes from the 
inner part of the disk in the Very High and High spec- 
tral states to the isotropically emitting ADAF in the 
Intermediate and Low states. 

We estimated the strength of irradiation required 
to reproduce the observed X-ray to optical flux ratio 
in A0620— 00, and found that we need an irradia- 
tion parameter (as defined by Eq. 4-2) of C ~ 0.004 
when the X-ray production is dominated by the disk, 
and C ~ 0.03 when the irradiating flux comes from 
the ADAF. This corresponds to a fixed disk height 
at the outer edge of H/R ~ 0.12, though this num- 
ber probably should not be taken at face value due 
to considerable uncertainty in the value of the X-ray 
albedo of the outer disk. However, for GRS 1124—68 
we estimated that the optical data taken during the 
VHS and HS requires an irradiation parameter a fac- 
tor of 3 smaller than in A0620— 00, corresponding 
to H/R ^ 0.07. We speculate that this difference be- 
tween otherwise nearly identical systems implies that 
outer disk fiaring is caused by warping. 
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Table 1 

Summary of the Binary Parameters for GRS 1124-68 and A 0620-00 



Source Name 



GRS 1124-68^ 



A0620-00 



Orbital Period, Porb (hr) 
Mass Function, /(M/Mq) 
Companion Type 
Black Hole Mass, m = M/Mq 
Mass Ratio, q = M^/M 
Inclination, i 

Outer Disk Radius} Rout (10^^ 
Distance, d (kpc) 



10.4 

3.01 ±0.15, 3.34 ±0.15'^ 
K3-5V 

6 n+i-^ 

"•U_l.O; 



cm 



0.133 zbO.Oigf 0.128 ±0.04' 
60l^o, 54^ 
1.9 ±0.16 
5.0 ± 1.1 



1+20° 
-15° 



7.75^ 

2.91 ±0.08^*2.72 ±0.06^ 
K5V 
10±7f 4.6± l.O"^ 
0.067 ± 0.01^ > 0.094^ 
37iP°g''^ 67±4°'^ 
1.2 ±0.27 
1.05 ± 0.4S''^ 



^Computed using equation (I 



Refer ences. — (a) From [EMN| and references therein, except where explicitly indicated ; (b) [Casares et al 
199 7f] (c) IShahbaz et al 1997t (d) Prosz et al. 1 99^; (e) pVIcChntock fc Remillard 1986|; (f) [Marsh et al. 1991 ; 
(g) IShahbaz et al. 1994| ; (g) paswdl et al. 1993| ; (k) [Shahbaz et al. 1999| ; (m) [Barret et al. 1996| . 



